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Abstract 


The  Multiple  Signal  Classification  (MUSIC)  algorithm  has  been  used 
widely  in  applications  with  spatially  distributed  sensor  arrays.  In  these 
applications,  the  algorithm  uses  the  phase  relations  of  the  signal  at 
different  sensor  pointing  angles.  Accordingly,  the  MUSIC  algorithm 
has  not  been  used  by  others  to  superresolve  target  positions  within 
the  main  beam  of  a  radar  with  a  single  narrow-beam  antenna  that  is 
step-scanned  in  angle  because  there  is  no  phase  shift  due  to  the 
scanning  of  the  beam.  We  have  developed  an  application  called  Scan 
MUSIC  (SMUSIC)  for  the  spatial  resolution  of  closely  spaced  targets 
with  a  stepped-scan  radar.  We  apply  SMUSIC  to  resolve  closely  spaced 
coherent  targets  when  a  one-dimensional  stepped  scanning  antenna  is 
used.  As  the  low-cost  enabling  radar  technology  (LCERT)  antenna 
points  at  different  angles  in  the  horizontal  (azimuth)  direction,  a 
step-scan  sensor  vector  is  formed.  In  this  report,  we  demonstrate  that 
SMUSIC  can  be  used  for  a  scanning  antenna  using  experimental 
LCERT  radar  data.  We  also  show  the  problems  encountered  when 
using  experimental  data  from  targets  exhibiting  constructive  and 
destructive  interferences. 
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1.  Introduction 


The  Rayleigh  criterion  limits  the  angular  resolution  of  an  antenna  aper¬ 
ture  of  size  D  to  the  value  9  =  A/D,  where  X  is  the  wavelength  corre¬ 
sponding  to  the  transmitted  frequency.  The  Multiple  Signal  Classification 
(MUSIC)  algorithm  [1],  which  uses  the  eigenvector  decomposition 
method,  is  a  superresolution  algorithm  widely  used  to  locate  closely 
spaced  multiple  emitters  (targets)  with  high  resolution  (smaller  than  the 
Rayleigh  criterion).  The  algorithm  is  often  employed  with  a  uniform 
linear  array  of  sensors  of  total  aperture  length  D  and  depends  on  the 
phase  progression  between  sensors  of  the  returned  signal.  In  contrast,  a 
single  antenna,  also  of  aperture  D  and  having  a  beamwidth  9,  which  is 
scanned  in  angle,  generates  no  phase  shifts  within  the  main  lobe  due  to 
the  target  direction.  Therefore,  the  MUSIC  algorithm  has  not  been  used 
by  others  to  superresolve  target  positions  within  the  main  beam  of  a 
single  narrow-beam  antenna.  In  this  report,  we  describe  an  application 
called  Scan  MUSIC  (SMUSIC)  [2]  that  is  used  to  improve  the  spatial 
resolution  of  closely  spaced  targets  with  a  stepped-scan  radar.  The  signal 
amplitude  vector  formed  by  the  response  of  the  antenna  as  it  is  stepped  in 
angle  replaces  the  vector  of  multiple  sensor  outputs  of  a  uniformly  linear 
spaced  array  of  sensors.  In  many  cases,  this  replacement  of  a  large  num¬ 
ber  of  sensors  by  a  single  rotatable  sensor  may  be  advantageous.  We 
demonstrate  that  superresolution  is  obtainable  in  this  stepped-scan  case, 
albeit  with  the  requirement  for  a  higher  signal-to-noise  ratio  (SNR)  than 
that  required  for  an  antenna  array  of  spaced  elements  having  the  same 
Rayleigh  resolution.  We  show  results  of  SMUSIC  using  the  data  obtained 
with  an  experimental  system  called  the  low-cost  enabling  radar  technol¬ 
ogy  (LCERT)  radar  [3]  as  the  radar  antenna  scanned  horizontally  in  the 
region  of  interest  (ROI). 
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2.  LCERT  Experimental  Data 

The  LCERT  instrumentation  radar  system  operates  at  Ku-band  (16  GHz) 
and  is  placed  on  a  truck  trailer  platform.  The  antenna  was  mounted  on  a 
pedestal  and  had  a  relatively  wide  elevation  beam  (cosecant-squared 
antenna)  and  a  narrow  azimuth  beamwidth  (less  than  1.5°).  The  radar  has 
two  operating  modes  (stationary  target  indicator  (STI)  and  moving  target 
indicator  (MTI))  and  different  modulation  modes  (linear  frequency 
modulation  (LFM)  and  phase  modulation  (PM)).  It  is  computer-controlled 
with  a  high  data  rate  recording  capability.  Table  1  shows  the  characteris¬ 
tics  of  LCERT  radar. 

The  experimental  data  for  this  SMUSIC  analysis  were  collected  at  an 
Aberdeen  Proving  Ground  (APG),  MD,  test  range  in  1997.  The  test  range 
is  a  flat,  grassy,  and  open  field.  The  coherent  returned  signals  were  col¬ 
lected  as  the  LCERT  antenna  scanned  horizontally  in  the  ROI,  which  had 
corner  reflectors  (targets)  within  the  main  beam  of  the  antenna.  Figure  1 
shows  the  antenna  used  in  the  data  collection. 

The  antenna  amplitude  pattern  was  obtained  using  an  aluminum  trihe¬ 
dral  (corner  reflector  or  target)  with  radar  cross  section  (RCS)  =  21.5  dBsm 
as  a  single  target.  This  was  placed  600  m  downrange  at  the  center  of  the 
radar  line-of-sight  axis.  The  antenna  amplitude  data  were  used  as  the 
array  manifold  vector  for  the  SMUSIC  algorithm. 

Two  other  aluminum  targets  (shown  in  fig.  2),  with  RCS  =  17.0  dBsm, 
were  placed  at  a  cross-range  position  600  m  downrange  from  the  radar. 

All  corner  reflectors  faced  the  radar.  The  antenna  continuously  scanned 
within  ±3°  of  the  ROI  with  respect  to  the  center  of  the  antenna  axis  (0°) 
(see  fig.  3).  The  targets  were  set  up  in  cross  range  as  shown  in  figure  4. 

The  experimental  scan  data  were  obtained  by  strobing  the  sampling  data 
from  the  slowly  scanning  antenna. 


Table  1.  LCERT  radar 
characteristics. 


Parameter 

Value 

Transmit  frequency  (GHz) 

16 

Transmit  pulsewidth  (ps) 

2.48 

Receive  pulsewidth  (ps) 

5.12 

Transmit  chirp  rate  (MHz/ps) 

8.065 

Transmit  start  frequency  (MHz) 

100 

Receive  start  frequency  (MHz) 

99.625 

Pulse  repetition  frequency  (kHz) 

8 

Analog-to-digital  samples  per  pulse 

256 

Pulses  per  coherent  processing  interval 

32 

Range  resolution  (m) 

10.8 

Number  of  coherent  processing  intervals 

1000 

Radar  mode  (STI  or  MTI) 

STI 

Modulation  (1-LFM,  2-PM) 

1 

Azimuth  angular  scan  swath  (°) 

8 

Scan  rate  (°/s) 

4 

2 


Figure  1.  LCERT 
antenna. 


Figure  2.  Corner 
reflectors  for 
experimental  setup. 


Figure  3.  LCERT 
experimental  setup. 


Figure  4.  Target 
placement  for  LCERT 
experimental  setup. 
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In  figure  4,  the  vertical  dotted  lines  are  the  actual  corner  reflector  posi¬ 
tions  at  600  m  downrange.  All  corner  reflectors  (targets)  are  within  the 
main  lobe  of  the  antenna  beam  since  the  cross  range  is  about  13.6  m  for 
the  antenna  beamwidth  1.3°.  All  the  runs  except  the  first  run,  calO.lcer, 
and  run  3  had  two  targets  placed  as  shown  in  the  figure  and  covered  as 
many  combinations  for  V2,  V3,  V4,  3/4  beamwidth  of  the  antenna.  For  runs 
10  to  13,  we  used  two  smaller  RCS  targets  (13.94  dBsm)  to  determine 
whether  SMUSIC  can  resolve  closely  spaced  targets  at  lower  RCS. 

The  radar  was  externally  calibrated  every  time  data  were  collected  by 
using  the  standard  method  of  checking  the  maximum  return  of  a  single 
corner  reflector  mounted  on  a  tripod. 

The  LCERT  radar  data  were  first  preprocessed  by  the  direct  sampling 
method  [4,3]  to  obtain  the  I  (in-phase)  and  Q  (quadrature  or  out-of-phase) 
components  as  well  as  the  magnitude  of  the  data,  which  are  suitable 
formats  for  the  SMUSIC  algorithm.  Figure  5  shows  the  layout  of  the  raw 
data  stored  in  a  data  file. 

The  data  contained  1000  coherent  processing  intervals  (CPIs).  Each  CPI 
consisted  of  32  pulses  (0, 1,  ... ,  31),  and  256  samples  for  each  pulse. 
Within  each  CPI,  the  frequency  was  chirped  (linear  frequency  modulated) 


Figure  5.  Raw  data 
format. 


CPI  1 


5 


with  a  bandwidth  of  20  MHz  for  each  pulse,  and  the  starting  frequencies 
for  the  pulses  were  set  at  increments  of  20  MHz  so  that  the  total  band¬ 
width  was  640  MHz.  Then  the  first  550  CPI  were  used  and  pulse  0  was 
selected  for  each  CPI  to  apply  the  SMUSIC  algorithm  (other  pulses  give 
the  same  results).  Figure  6  is  a  block  diagram  of  the  preprocessed  LCERT 
data  before  the  SMUSIC  process. 

In  figure  7,  data  from  the  first  run,  calO.lcer,  show  the  antenna  pattern 
obtained  when  the  LCERT  radar  antenna  scanned  within  ±3°  of  the  ROI 
for  a  single  reflector.  Figures  8  and  9  are  examples  of  the  responses  that 
exhibit  two  types  of  interference  from  constructive  and  destructive  inter¬ 
ference.  Figure  8  shows  the  target  response  of  two  reflectors  within  the 
main  lobe  of  the  antenna  that  constructively  interfere.  Figure  9  shows  the 
target  response  of  the  same  two  reflectors  when  they  destructively 
interfere. 


Figure  6.  Block 
diagram  of 
preprocessing  of 
LCERT  raw  data. 
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Figure  7.  LCERT 
antenna  pattern. 
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Figure  8.  Constructive 
target  response. 
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Figure  9.  Destructive 
target  response. 
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3.  SMUSIC  Process 


Figure  10  depicts  the  processing  sequence  of  the  SMUSIC  algorithm.  In 
this  figure,  the  LCERT  block  provides  the  LCERT  data  collected  at  APG. 
The  data  samples  are  the  data  centered  with  respect  to  the  center  of  the 
antenna  pattern.  The  data  were  then  truncated  between  +3°;  these  are 
called  truncated  data.  Finally,  the  data  samples  were  taken  from  every  6th 
sample  of  the  truncated  data  between  ±1.5°  to  be  used  for  the  SMUSIC 
computation.  These  steps  were  done  in  the  PRE-format  block. 

In  general,  we  can  write  the  correlation  matrix  of  a  complex  random 
vector,  x,  as  R  =  e[xxH],  where  x  is  the  vector  of  observed  data  and  is 
referred  to  as  a  snapshot,  and  the  elements  of  a  snapshot  are  called 
samples.  e[  ]  is  the  expected  value  or  ensemble  average  operator,  and  H  is 
Hermitian  transpose.  Then  the  data  vector,  x,  with  N  samples,  was  used 
to  estimate  the  correlation  matrix.  Since  there  is  only  a  single  snapshot 
(scan)  and  the  target  response  is  coherent,  Rxx  will  be  rank  deficient.  In 
fact,  the  rank  of  Rxx  is  equal  to  1  (a  proof  of  this  property  is  given  by 
Shan  et  al  [5] ).  The  following  simple  example  demonstrates  this  problem. 

Let  a  vector  x  =  [1  3  5  7  4]^.  Then 
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Then  the  rank  of  Rxx  is  equal  to  1  because  the  columns  of  Rxx  are  linearly 
dependent.  To  avoid  the  rank  deficient  problem,  Evans  [6]  developed  a 
technique  called  forward /backward  spatial  smoothing.  For  more  infor¬ 
mation  and  a  further  explanation  of  this  technique,  see  Van  Trees  [7].  In 
this  application,  we  only  use  forward  spatial  smoothing  to  obtain  the  full 
rank  of  Rxx  and  also  decorrelate  the  coherent  signals.  We  call  this  the 
forward  subvector  averaging  technique  and  use  it  to  estimate  the  correla¬ 
tion  matrix  as  described  below.  The  nonsingularity  of  the  sample  correla¬ 
tion  matrix,  Rxx ,  is  crucial  to  an  application  of  the  eigenvector  decompo¬ 
sition  method;  the  technique  guarantees  the  nonsingularity  of  the  sample 
correlation  matrix  of  the  signals. 


Figure  10.  Block 
diagram  of  SMUSIC 
process. 
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Let  x  be  an  N  x  1  input  data  vector.  Define  subvectors  xw  for  1  <  i  <  p 
(number  of  subvectors),  to  be  a  vector  containing  Q  consecutive  samples 
from  x  with  Q  <N,  starting  at  sample  i  as  shown  in  equation  (1): 


x(D  = 
x(2)  = 


xl'  x2'  ,  XQ 


x2'x3'  'XQ  +1 


(1) 


r(p)  = 


xp,  Xp  +  p 


xp  +  Q- 1 


Since  there  are  N  samples  in  x,  the  possible  values  of  p  and  Q  are  related 
to  Nby 

p  +  Q-l=N.  (2) 

Next,  let  the  subvector  correlation  matrix  be  of  x®  x^H.  Finally, 
define  the  estimate  correlation  matrix  Rxx  to  be  the  average  of  p 
subvector  correlation  matrices  as 

=  R8=1.L«*«"  ,  (3) 

computed  in  the  Rxx  block.  Note  that  the  matrices  R^  ^  and  Rxx  are  of  the 
order  of  Q  x  Q ,  while  R  is  of  the  order  of  N  x  N .  So,  by  the  subvector 
averaging  technique,  the  dimension  of  the  sample  correlation  matrix  is 
effectively  reduced.  The  following  example  demonstrates  this  technique. 

Let  x  be  the  input  data  vector  as  used  above,  let  p  =  3,  so  Q  =3;  then 

x<1)=  [13  5]7,x(2)  =  [3  5  7]7,andx(3)  =  [5  7  4]r  . 

Then  Rxx  is  computed  using  equation  (3)  as 


R 


xx 


11.6667  17.6667  15.3333 

17.6667  27.6667  26.0000 

15.3333  26.0000  30.0000 


and  the  rank  of  Rxx  is  3,  which  is  full  rank.  We  see  that  the  original 
correlation  matrix  has  the  dimension  5x5,  and  the  estimate  correlation 
matrix  using  the  subvector  averaging  technique  has  the  dimension  3x3. 
The  trade-off  with  this  technique  clearly  shows  that  the  dimension  is 
effectively  reduced  but  the  rank  is  increased. 


9 


If  the  estimated  correlation  matrix  is  nonsingular,  the  distinct  eigenvalues 
and  corresponding  eigenvector  space  from  the  estimate  Rxx  are 

>  ^2 , Aq}  and  { Vy  V2, . . .  Vq}  .  (4) 

The  computation  is  carried  in  the  EIVV  (eigenvectors  and  eigenvalues) 
block  (see  fig.  10). 

From  the  argument  above,  if  there  are  M  targets,  the  rank  of  Rxx  will  be 
M.  It  follows  that  the  minimal  eigenvalues  of  Rxx  are  j  A;V1  + 1  >  +  2, . . .  , 

>  Aq},  and  the  eigenvectors  corresponding  to  the  minimal  eigenvalues  are 
{ vM  +  y  vM  +  2,  ■  ■  ■  t  Vq}.  These  eigenvectors  span  the  noise  subspace.  The 
orthogonal  complement  of  the  noise  subspace  is  referred  to  as  the  signal 
subspace.  The  signal  subspace  denotes  Us  ,  and  the  noise  subspace 
denotes  UN .  The  dimension  of  Us  is  Q  x  M  and  UN  is  Q  x  (Q  -  M).  In 
this  experiment,  since  we  know  the  number  of  targets  in  the  field,  M  =  2, 
we  use  this  information  to  separate  the  eigenvector  space  into  a  signal 
eigenvector  with  dimension  Q  x  2  and  noise  eigenvector  with  dimension 
Q  x  (Q  -  2)  subspaces  in  the  Signal/noise  subspaces  block. 

We  compute  the  SMUSIC  spectrum  using  equation  (5)  in  the  Compute 
inverse  distance  block.  The  SMUSIC  algorithm  uses  the  orthogonality 
property  that  exists  between  the  signal  eigenvector  subspace  and  the 
noise  eigenvector  subspace  from  the  sample  correlation  matrix  to  estimate 
the  position  of  each  target  center.  This  can  be  accomplished  by  computing 
the  inverse  of  the  distance  of  a  searching  or  true  direction  or  an  array 
manifold  vector  at  angle  9  and  the  noise  eigenvector  subspace.  The  peaks, 
computed  in  the  Peaks  of  SMUSIC  block,  from  the  SMUSIC  spectrum  are 
the  estimated  target  locations 


W0)=(^)HUNUNHv(0)j  ,  (5) 

where  v(9),  that  is,  the  scanning  antenna  pattern,  is  the  steering  vector  or 
array  manifold  vector  for  all  values  of  9. 

In  the  SMUSIC  computation,  we  used  the  LCERT  antenna  pattern  as  the 
array  manifold  vector  to  compute  the  spectral  SMUSIC.  The  array  mani¬ 
fold  vector  for  this  scanning  antenna  in  this  report  is  a  magnitude  two- 
way  scanning  antenna  pattern.  Unlike  the  ordinary  MUSIC,  which  mostly 
uses  a  uniformly  linear  array  manifold  vector  for  the  computation,  the 
steering  vector  of  the  SMUSIC  does  not  depend  on  the  phase  relationship 
from  one  step  scan  to  another  step  scan. 
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4.  Results 


In  this  section,  we  show  the  LCERT  experimental  data  collected  at  APG  in 
section  4.1  and  the  results  of  SMUSIC  in  section  4.2. 


4.1  LCERT  Data 

The  magnitudes  of  the  target  responses  in  decibels  versus  azimuth  scan 
angle  in  degrees  for  all  runs  are  shown  in  figures  A-l  to  A-13  in  appendix 
A.  Figure  A-l  is  the  LCERT  antenna  pattern  used  as  the  array  manifold 
vector  for  the  SMUSIC  computation.  From  this  figure,  we  see  that  the 
noise  from  the  response,  including  the  system  noise  and  the  background 
noise,  is  small  compared  to  the  signal  response.  From  figures  A-2  to  A-13, 
we  see  that  the  responses  exhibit  varying  types  of  interferences  from 
constructive  to  partially  constructive  to  destructive  interferences.  Figures 
A-2,  A-4,  A-6,  A-8,  and  A-13  show  the  target  responses  of  two  reflectors 
within  the  main  lobe  of  the  antenna  that  constructively  interfere.  Figures 
A-3,  A-5,  A -7,  and  A-9  to  A-12  show  the  targets  responses  of  the  same  two 
reflectors  when  they  fully  and  partially  destructively  interfere.  These 
types  of  interferences  affect  the  performance  of  SMUSIC. 


4.2  SMUSIC  Results 

LCERT  data  included  two  different  cases:  constructive  and  destructive 
interferences.  For  the  constructive  interference,  SMUSIC  resolved  the 
closely  spaced  targets  when  the  targets  were  placed  apart  at  distances  less 
than  the  beamwidth  of  the  antenna.  On  the  other  hand,  for  the  destructive 
interference,  SMUSIC  was  not  able  to  resolve  the  closely  spaced  targets 
within  the  beamwidth  of  the  antenna.  In  attempts  to  resolve  targets 
having  destructive  interference,  we  use  the  phase  information  of  the  data 
as  well  as  the  phase  information  of  the  manifold  vector.  We  used  complex 
arithmetic  for  the  SMUSIC  process.  However,  SMUSIC  was  still  unable  to 
resolve  the  targets. 

Figures  B-l  to  B-12  in  appendix  B  show  the  results  of  SMUSIC  for  LCERT 
data  runs.  In  this  SMUSIC  application,  we  only  used  the  magnitude 
information  of  the  scanning  antenna  and  the  target  response  for 
superresolving  the  targets  within  the  main  lobe  of  the  antenna  beam. 
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5.  Conclusions  and  Future  Investigations 

We  have  demonstrated  that  a  superresolution  algorithm  like  MUSIC  can 
be  extended  for  a  scanning  antenna  even  though  there  is  ideally  no  phase 
variation  for  the  scanning  antenna  within  the  main  beam  of  an  antenna. 
We  have  developed  an  application  called  SMUSIC  to  apply  to  the  scan¬ 
ning  antenna  that  used  the  magnitude  of  the  target  response  to 
superresolve  two  closely  spaced  targets  within  the  main  beam  of  the 
LCERT  radar  antenna.  SMUSIC  resolved  the  closely  spaced  targets  when 
they  interfered  constructively;  i.ev  they  were  within  30°  or  less  phase  shift 
of  each  other. 

We  showed  the  results  of  SMUSIC  for  all  the  cases  of  partially  or  fully 
constructive  interference.  We  conclude  that  SMUSIC  can  be  used  for  real 
data  with  a  narrow-band  signal  and  that  it  performs  beyond  the  limit  of 
the  Rayleigh  criterion.  It  can  resolve  the  target  locations  in  which  the 
target  separation  is  less  than  one-third  of  an  antenna  beamwidth. 

For  the  fully  and  partially  destructive  interference  cases,  SMUSIC  may 
not  be  able  to  resolve  the  targets  within  the  antenna  beamwidth.  This 
problem  needs  further  investigation. 

An  advantage  of  a  scanning  antenna  is  that  it  avoids  the  requirement  of  a 
phase  calibration  of  a  linear  array  of  sensors.  A  single  sensor  can  be 
cheaper  than  making  multiple  sensors  with  the  same  aperture  length. 

A  more  comprehensive  Monte  Carlo  study  is  needed  to  analyze  the 
performance  of  the  SMUSIC  algorithm  for  a  scanning  antenna  in  terms  of 
the  antenna  beamwidth,  mean-square-error  (MSE),  probability  of  success¬ 
ful  resolution  (PSR)  in  a  function  of  the  number  of  step  scans,  numbers  of 
subvector  averaging,  and  SNR  for  two  coherent  targets  with  different 
degrees  of  correlation. 
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Appendix  A.  LCERT  Magnitude  Response 
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Appendix  A 


Figure  A-l.  LCERT 
antenna  pattern  (skip 
data  point  (SDP)  =  6). 


Figure  A-2.  LCERT 
data  run  1  (SDP  =  6) 
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Appendix  A 


Figure  A-3.  LCERT 
data  run  2  (SDP  =  6). 
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Figure  A-4.  LCERT 
data  run  4  (SDP  =  6). 
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Appendix  A 


Figure  A-5.  LCERT 
data  run  5  (SDP  =  6). 


Figure  A-6.  LCERT 
data  run  6  (SDP  =  6). 
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Appendix  A 


Figure  A-7.  LCERT 
data  run  7  (SDP  =  6). 


Figure  A-8.  LCERT 
data  run  8  (SDP  =  6). 
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Appendix  A 


Figure  A-9.  LCERT 
data  run  9  (SDP  =  6). 


Figure  A-10.  LCERT  20 

data  run  10  (SDP  =  6). 
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Appendix  A 


Figure  A-ll.  LCERT 
data  run  11  (SDP  =  6). 


Figure  A-12.  LCERT 
data  run  12  (SDP  =  6). 
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Appendix  A 


Figure  A-13.  LCERT 
data  run  13  (SDP  =  6). 
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Appendix  B.  SMUSIC  Outputs 


25 


Appendix  B 


Figure  B-3.  SMUSIC 
output  for  run  4 
(SDP  =  6,  SVA  =  4). 


Figure  B-4.  SMUSIC 
output  for  run  5 
(SDP  =  6,  SVA  =  4). 


Appendix  B 


Figure  B-7.  SMUSIC 
output  for  run  8 
(SDP  =  6,  SVA  =  4). 


Figure  B-8.  SMUSIC 
output  for  run  9 
(SDP  =  6,  SVA  =  4). 
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Appendix  B 


Figure  B-ll.  SMUSIC 
output  for  run  12 
(SDP  =  6,  SVA  =  8). 


Figure  B-12.  SMUSIC 
output  for  run  13 
(SDP  =  6,  SVA  =  8). 
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